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O ' Abstract 

VO . 

■ The behavior near the singularity of an isotropic, homogeneous 

O cosmological model with a viscous fluid source is investigated. This 

turns out to be a relaxation dominated regime. Full extended irre- 
^Jy versible thermodynamics is used, and comparison with results of the 

truncated theory is made. New singular behaviors are found and it is 
. £h \ shown that a relaxation dominated inflationary epoch may exist for 

^ ' fluids with small heat capacity. 
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There are a many processes capable of producing important dissipative 
stresses in the early universe. This includes interactions between matter and 
radiation quarks and gluons 0, different components of dark matter [f|], 
and those mediated by massive particles. It happens also due to the decay 
of massive superstrings modes ||, gravitational particle production || || 
and phase transitions. Phenomenologically, these processes may be modeled 
in terms of a classical bulk viscosity. The dynamics of the viscous pressure 
term is ruled by a transport equation, and we use the expression given by 
the Extended Irreversible Thermodynamics theory @ [g]. 

a + T& = -^ H -i Ta ^ H + ^-^-^\ (1) 

where we are assuming a spatially flat Robertson- Walker metric, a is the bulk 
viscous pressure, H = a/a the Hubble variable, r is the relaxation time, ( is 
the bulk viscous coefficient, T is the temperature, and e = corresponds to 
the truncated theory while e = 1 corresponds to the full theory. Equation ([I]) 
yields a causal and stable behavior due to the memory mechanism provided 
by the relaxation time. We assume that £ = ap m , r = (/p and T = np r , 
where a > 0, m, n > and r > are constants [f|. 

We will examine in this Letter the evolution of the dissipative universe 
near a singularity, which is a relaxation dominated era. In effect, the Einstein 
equations are 

H 2 = \ P (2) 

H + ZH 2 = l -{p- V -a) (3) 

where p is the energy density, p is the equilibrium pressure and we take the 
equation of state p = (7 — l)p with a constant adiabatic index < 7 < 2. 
For a singularity in to, we will consider two types of leading behaviors when 
t — > to, H ~ 1/At n , where n > is a constant, and H ~ l/(AtlnAt). For 
H ~ 1/Ar, dociJ 2 ~ l/(At) 2n for n > 1, and a oc H ~ l/(At) n+1 for 
n < 1. So, |<x/(C#)| ~ Af^ 2 ™- 1 ) for n > 1, and |cr/(C#)| ~ At 2 ™ 71 " 1 for 
n < 1. Also f/r ~ ~ T/T ~ 1/At. Thus we find that the first term 
of the left hand side of (0) becomes much smaller than the rest of the terms 
provided that m > 1/2 and r 3> |A£|. Similar considerations are valid for the 
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second type of leading behavior. Then, using ([[]), (0) and fl|), and neglecting 
the viscous term a, we obtain 



H-e(l+r)— + 3{ 7 + ^[l-(l+r) 7 ]}M+ i (e 7 -2)if 3 = (4) 

We investigate the solutions of this equation with one of these two types 
of behaviour near the singularity, as the solutions of the non-approximated 
equation will have the same leading behavior for At — > 0. The study of 
these solutions for large time will be made elsewhere (see also [7]). As the 
relaxation regime has already been considered in the truncated theory []T0f| , 
we will concentrate on the case e—1. We also note that an equation with the 
same structure arises when considering the semiclassical Einstein equation 
with vacuum polarization terms near the singularity |Q|. 

To verify the consistency of our approximation, we check first whether 
equation (|J) has solutions of the form a = aoAt u . In effect there is a region 
of the parameter plane (7, r), such that two families of solutions of this form 
exist, with v given by 



v± = W^T) {" [1 + (1 " r)7] ± t (1 + (1 " r)7)2 ~ 4(2 " 7)(r " 1} ' 

(5) 

This region is determined by the requirement D > 0, where D is the dis- 
criminant in (|^), and this condition can be satisfied for 7 < 2. On the curve 
D = there is only one family of solutions, and for 7 = 2 there is also one 
family of solutions with v-i = 2(1 — r)/(3(3 — 2r)). 

To investigate further the two-parameter families of solutions of it 
is convenient to make the change of variable H = y -1 ^, which turns this 
equation into 

y + -[l + (l- r) 7 ] y-^ij + ^(2 - ^ry^r = (6) 

In the case that 1 + (1 — r) 7 7^ 0, r ^ 1 we solve (|(|) applying a generalization 
of the transformation used in 11121 



1 3 1 

-y~ , dr\ — — [1 + (1 — r) 7 ] y~~dt 



r-V ' 2 
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which linearizes (|6]) 

d 2 z ^ dz ^ _ q 
drf dr) 

where j3 = (2— j)(r— 1) / [1 + (1 — r)^} 2 . Thus, we obtain the general solution 
of (0) in parametric form 



r — 1 



H(rj) 



A ^- 3[l + (l-r) 7] y^ (8) 

This parametric form is very useful to obtain the singular behavior in the 
region D < 0, where no explicit solution is available. In this region r > 1, 
so that H diverges when z — > 0. Moreover, as this region corresponds to 
the oscillating or critically damped regime of the oscillator described by (|7|), 
when 7 < 2 z(rj) oc Arj for Arj —>■ 0. Then we find H ~ C / 'At" or 

a(t) ~ a exp (KAt 1 ^) , At -> (9) 



where ao and X are arbitrary integration constants. This behavior exhibits 
a singularity at a finite value of the scale factor. 

In the case that 1 + (1 — r)7 = and 7 < 2, equation (^) describes the 
motion of a point particle in a potential V(y) oc y 2 ( r ~ 1 )/ r and reduces to 
quadratures 

At = = / dH m (10) 

where E is an integration constant and G = 9(2 — 7 )(1 + 7 ) 2 /(8 7 ). This 
solution may be expressed in terms of the hypergeometric function, and we 
find that the singular behavior is also (|9|), where r = (1 + 7 )/ 7 . 

On the other hand, for D > 0, explicit two-parameter families of solutions 
of (|G|) exist provided that the constrain 



< 2 "* -'2- 1 -) (11) 



[1 + (1 - r) 7 ]' 



r 
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is satisfied. They yield a representative sample of two-parameter families of 
solutions in this region, and have the form 



a(t) = a 



\At\ 



+ K 



r ^ 1 



(12) 



a(t) = a | In \At\ + K\ 



(13) 



where u y = 2(1 — 2r)/[3(l + (1 — r)j)]. We find in this way four kinds of 
singular behaviors, as shown in the following table: 



a(t) 


lim a(t) 


Parameter range 


At v + 





< r < 1 and r > 3/2 


At u - 


oo 


< r < 1 and 1 < r < 3/2 


llnAtl"- 





r = 1 


exp(iTAt 2 v i ) 


a 


r > 1 



Finally, for 7 = 2, equation (0) also reduces easily to quadratures 



where S = 3r(3 — 2r)/(2r — 2), and this solution be expressed in terms of 
the hypergeometric function. We find that their singular behavior is At V2 for 
< r < 1, At~ 2/3 for r = 1, and (|) for r > 1. 

We have investigated the singular behavior of the full causal model for 
m > 1/2. Similarly to the truncated causal model, we find that it is domi- 
nated by the relaxation terms of the transport equation. However we have 
found that a greater variety of singular behaviors appear in the full model, 
and it is the dependence of the energy density on the temperature which de- 
termines the kind of singular behavior of two-parameter families of solutions. 

Big-bang and explosive singularities already appeared in the truncated 
theory. Following the same steps as before and assuming H = vj At we find 



v± = (1/3) [ 7 ±(7 2 + 4) 
for < r < 1, while 1/4. 



1/2' 



> 



(cf. 0). We note that < u+ < v + < 2/ (37) 
_ 4/3 > u + and u + > 2/(37) for r > 3/2. Thus 
we can state that for a fluid with large heat capacity, the rate of expansion 
of the scale factor after a Big-bang singularity is smaller than calculated in 
the truncated theory which is turn smaller than the perfect fluid rate. In 
this case, there are particle horizons. However, when the heat capacity is 
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small, the rate of expansion turns out much larger for a viscous fluid than a 
perfect fluid. Moreover, it turns out to be a relaxation dominated inflationary 
regime. 

On the other hand, z/_ < £/_ < for < r < r c and P_ < za_ < for 
r c < r < 3/2, where r c depends on 7. Thus, explosive singularities may be 
stronger or milder depending on the heat capacity being large or small. A 
finite scale factor singularity arises for r > 1, and a logarithmic explosive 
singularity occurs for r = 1. Neither of them appear in the truncated model, 
though finite scale singularities has been found in noncausal models for m > 
1/2 0. 
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